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A theoretical view of practical problems in interferometry 
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Interferometry is discussed in terms of the representation of the source. In particular, scale- 
invariant 1-d hydrodynamics is revisited, and extended to the case of unequal transverse masses. 
It is argued that that kaon emission occurs over a short time interval. Exact results for models of 
two- and three-dimensional flow are presented, which exhibit altered scaling laws. Such qualitative 
trends, together with other observables, are vital if one is to draw conclusions about the source M. 



I. INTRODUCTION 



On 
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The discovery a quark-gluon plasma will to some extent be linked to the measurement of the geometric size and 
other macroscopic characteristics of the reaction zone. An important tool in accomplishing this is interferometry. 
As proposed by Hanbury-Brown and Twiss to measure stellar sizes, and as the GGLP effect in (low energy) nuclear 
physics, it provides for a straightforward determination of the size of the "hot source." In high-energy collisions, on 
the other hand, the correspondence between measured quantities and parameters of the emitting system is less clear 
H . For example, correlations on the same scale as the size of the emitting system make an interferometric study of 
■ the correlations inside the matter and a measurement of the external geometry exceedingly difficult . 

For A-A collisions, even in the simplest scenario of hydrodynamic evolution, the relation between the inclusive 
H ■ one- and two-particle spectra and the parameters of the emitting system does not follow the classical scheme of 
^ ' interferometry Since no direct inverse solution exists, and the question of what is measured by the two-pion 
correlation function becomes nontrivial. The problem for interferometry is thus to determine the parameters of a 
judiciously chosen density matrix p for the system, by measuring the inclusive cross-sections dNi/dk, dN%/ dkidk2, .... 
Here, we shall restrict ourselves to particle sources that may be described in terms of (semi-) classical one-particle 
distributions, or even by only a few macroscopic parameters. Practically, current dynamical models produce precisely 
\Q , such sources. However, even with these restrictions, the two-particle spectra are not easily decoded. We therefore 
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suggest in Sect gl| an illustrative "step-by-step" strategy: Essentially, the data are first confronted with the simplest 
physically motivated model, and discrepancies are then used show directions for improvement. Ideally, this requires 
that the model should have an "adequate" (semi-) analytical solution for the correlator, allowing one to obtain the 
model parameters precisely and directly, without resorting to intermediate fits to the data. 

Recent SPS experiments |||| have observed m± scaling of longitudinal radii, as well as the y-dependence of the 
correlator width predicted by boost-invariant 1-d hydrodynamic expansion 0. This provides a strong reason for 
taking macroscopic collective behavior (i.e., hydrodynamics) seriously and using interferometry to study it - at least 
at SPS energies. Here, we shall therefore extend the calculation of Ref. to give the (analytic) solution in the case of 
unequal transverse momenta, and illustrate why the data support a relatively rapid freeze-out picture for both pions 
and kaons. We also discuss natural extensions of the model to one that might be more appropriate at the AGS: (1) 
1-d Landau hydrodynamics, and (2) some exactly calculable examples with two- and three-dimensional flow. In these 
cases, the scaling law is altered, and we argue generally that such qualitative trends, together with other observables, 
are important in establishing an adequate model of of the source, and only then the value of its parameters. 

Before turning our attention to these matters, however, we discuss a related issue: there is no full consensus on 
what mathematical formulation is adequate for applications of interferometry. For example, an analysis of polariza- 
tion effects in photon interference M indicates that the languages of locally defined states and Wigner phase-space 
distributions lead to different answers; there is a practical need to resolve such controversies. We shall therefore begin 
in Sec. [n] by briefly formulating the problem of interferometry for an expanding source at the operator level. This 
emphasizes the quantum nature of the problem, and allows us to solve (correctly) the basic problem of one- and 
two-particle propagation without having to address the unknown nature of the source. With these expressions in 
hand, we discuss three examples of source input. 
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II. PRELIMINARIES 



Let |in) be any one of the initial states of a system emitting a pion field <p(x) at time t c . This field is detected later 
by an analyzer, tuned to the measurement of the momentum of the free pion; its eigenfunctions are thus plane waves 
ft{x) with x° > t c , and the corresponding annihilation operators for momentum k in the final state are given by 



d 3 xfUx)id x <p{x) 



(2.1) 



where (ad x b) = a(d® b) — (d® a)b. The operator A^ describes the effect of a detector far from the point of emission, 
so, by definition, the pion is detected on-mass-shell. The inclusive amplitudes to find one pion with momentum fc, 
or two with momenta fci, k 2 in the final state are (X\ AtS \in) and (X| Ag A^ S^in), respectively. Here, S is the 
evolution operator after freeze-out, and the \X) form a complete set of all possible secondaries. One obtains the one- 
and two-particle inclusive spectra by summing the squared modulus of these amplitudes over all (undetected) states 
\X), and averaging over the initial ensemble p. 
An analyzer tuned to k for singles, or fci, k 2 for pairs, is represented by the operators 



N(k) = A\A % , and NihM = A^A^A^A^ = - <5(fci - fc 2 )) , 

respectively. In the simplest description, the pion field ip(x) reaches the detector after free propagation 

<p(x) = / effi M (y) G rct (a; - y) d v <p{y) , 



(2.2) 



(2.3) 



where the field <p(y) and the space of states in which the density matrix acts are defined on the freeze-out hypersurface 
S c . The pion Fock operator (2.1) may now be expressed in terms of the initial fields, and simplified using the explicit 
form of the free pion propagators: 

&A S S = 9(x° - j, ) J dS^y) %(y)i T y viv) ■ (2.4) 

Clearly, the analyzers perform an on-shell Fourier expansion of the initial data. Collecting these results, we may 
rewrite the number operators for single and pairs as 



N(k) = 8(x"-y») dZMdZ u (y 2 ) 



*-> f , 
fk(yi)id yi ^{y 



0(Va)id V2 /|(j/: 



N(k u k 2 ) = d(x v -y u ) / dS^y 1 )dS v (y 3 )dE p {y4.)dSx(y2) 

> f_L i — y V 



d yi d y2 TMft (y: 



'fci 



(2.5) 



(2.6) 



respectively. The field operator products in these equations are now to be averaged over the density matrix of initial 
states. The interferometry problem requires that the basis the final quantum states of the preceding stage of evolution 
should be described explicitly. We give three examples: 

1. Emission from two one— dimensional cavities: We have two sets of the eigenstates 



for Ln — a < y < Lpj + a , where £jv = ±L 



= , otherwise. (2-7) 
The field decomposition ip(y) — N a p ,N 4>p,N(y) acquires an additional index enumerating the cavities, and 



states belonging to different cavities are independent: [a p> jsr, a p , N ,} 
distribution for when the walls are removed at yo — is obtained as 



S pp 'Snn'- 



N p 



(ujk + ujp) 2 sin 2 (p — k)a 
AiOkCOp na(p — k) 2 



From Eq. (2J5), the one-particle 



(2. 



where n(p, N) — (aj n o, p .n) is the boson occupation number in cavity N. 
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For the two-particle spectrum there is interference between four amplitudes: two when each particle is emitted from 
a dif fere nt cavity, and two more when both particles originate from the same cavity. For the statistical average of 
Eq. (2J3) we choose 



'&{ 0.3040,2} = n(l)n(3) 



^PlP2^A r iA r 2^P3P4^3A r 4 + < ^PlP4^A r iA r 4^P3P2^ 3 JV 2 



where we have introduced the notation i = (pi, Ni). Then 



(N( kl , fc 2 )> = (N(h)) (N(k 2 )) + V L+g^ +c pi )(^ 2 +„,,) 



1,3 



(2.9) 



(2.10) 



2^pi^p 3 



n(l)n(3) 



smjpx - fci)q sin(pi - fc 2 )o sin(p 3 - fci)a sin(p 3 - fc 2 )Q c -i(L 1 -L 3 )(fc 1 -fc 2 ) 



(Pi-fci) (pi-fo) (P3-fci) {P3-k 2 ) 



where Lj = L(Nj) — ±L, and momentum originates from the cavity iVj. The two terms in the sum with JVx ^ -/V3 
lead to the usual interference term, with Afc ~ 1/(2L). The average of the product of Fock operators (2.9), i.e., 
correlations in the emitting system, depends only on "internal" variables which identify the stationary states in 
the cavities before they are opened. On the other hand, the sin-functions serve as the projectors of these states 
onto the states of free propagation in which the particles are detected, and require additional physical input. First, 
detector resolution should not allow one to determine which cavity the particle originated from, i.e., we do not have 
fci — k 2 \ L 3> 1. This condition corresponds to the Rayleigh criterium: we are not able to construct an "optical" 
image of source. Secondly, in most physical situations L > a, so that \k% — k 2 \ a <C 1. This inequality removes the 
momentum integral in Eq. (|]^), so that a measurement of the one-particle spectrum does not allow for a determination 
of the cavity size. 

2. Emission from many cavities and hydrodynamics: To generalize the two-cavity model to the freeze-out of 
an extended hydrodynamic system, we consider a continuous set of decaying thermal cells, defined on a freeze-out 
surface T(x) = T c and Doppler shifted to account for cell motion 0. As in the previous example, the initial quantum 
states in one cell are assumed to commute with states in a different cell. This additional information implies the 
existence of two length scales, and allows us to use a quasi-classical picture. Within a cell, the system is governed 
by the short-range dynamics conveniently described in momentum language and if the interactions are switched 
off, the particles begin to propagate freely with the frozen momentum distribution defined within the cell. To the 
lowest approximation of the dynamical interactions and in agreement with the freeze-out concept, these particles are 
on-mass-shell. In terms of an auxiliary "emission function" J(k\,k 2 ), 



JT,,. 2 



one obtains for the one- and two-particle spectra Q 

,0^1 n.n dN; 



dk 



J{k, k) , and klk\ 



dk\dk 2 



J{ki, ki)J(k 2 ,k 2 ) +Re 



J(k 1 ,k 2 )J{k 2 ,k 1 ) 



(2.11) 



(2.12) 



respectively. The latter equation is explicitly symmetric in xi and x 2 , as well as in fc x and k 2 . It reveals interference 
of the two amplitudes for the two-pion state (consisting of the two on-mass-shell pions with quantum numbers k\ 
and fc 2 ) to be emitted by the two sources. These sources are located at points x\ and x 2 and have true thermal 
distributions, n(fcj • u{xi)), with on-shell particles. 

3. Wigner formalism: Simulations of heavy ion collisions typically produce a final state in terms of Wigner 
distributions, so it is of practical value to examine the two-particle spectrum in terms of them. Instead of Eq. ([2.9]), 
we may express the independence of the dynamical processes in the emitting system via 



(2.13) 



which implies that the distance z = x — y between the points in one correlator does not exceed a correlation length a cor 
which, in turn, is much less than the full size L of the system. In the first term, for example, x,y ~ R%, x — y < a cor 
and x' , y' ~ R 2 , x' — y' < a cor , while R\ — R 2 3> a cor . Under these conditions, it is possible to think of the coordinates 
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Rj as labels of domains of non-vanishing correlation, while the Fourier transform over z yields the local spectrum. 
Thus, changing variables^ x = R + z/2, y = R— z/2, and introducing the Wigner function Af 



^{R+ Z -)^R- Z -))= / d 3 sM(R,s)e^ 



Eqs. (fOp, Q2.13D , and ( |2.14| ) yield 

{N(kiM)) = j dRidz 1 ds 1 N{RiJi)e l{si - kl)zi J dR 2 dz 2 ds 2 Af(R 2 , s 2 )e l{s2 ^ k2)z2 + 



(2.14) 



(2.15) 



dRi dR 2 dz\ dz 2 ds\ ds 2 J\f(Ri, si)Af(R 2 , s 2 )e' 



If the range of integration a cor (which, strictly speaking, depends on R) for the z-variable is sufficiently large compared 
to 1/fci, 1/&2, the integrals over z\ and z 2 approximate to delta-functions, and the common expression for the correlator 
C 2 in terms of functions N{[k\ + k 2 )/2, R) results [||. Physically, however, since the (p{x) are free, there are no states 
with s 2 7^ to 2 . We thus face the problem of understanding the origin of a thermal distribution of off-shell pions, 
instantly frozen, and detected later in physical on-shell states. 2 The variables of integration do not correctly establish 
the correspondence between initial and final states. Moreover, the replacement W((fci + k 2 )/2,R) undermines the 
concept of independent emission of the two sources - a necessary condition for interference. Each source "knows" 
about the momenta of both pions, and one may reasonably suspect that emission from spatially separated volumes is 
not independent. 



III. INTERFEROMETRY AND SOURCE MODELS 



1. One— dimensional flow: We revisit interferometry for scale invariant 1-d hydrodynamic motion; evidence for 
such sources at SPS J5B comes from agreement with the theoretically predicted m t - and y-dependence of the 
"visible longitudinal size:" 



Ri 



coshy V ffh 



(3.1) 



where y is the pair mid-rapidity, and to 2 = m 2 + p 2 . Physically, the correlator measures the effective size of the fluid 
slice which forms the observable spectrum at a given rapidity. The larger the transverse momentum of a particle in 
the fluid, the more that particle is frozen into collective longitudinal motion: the subset of particles with high p± will 
have a very narrow distribution over k tl . 

The dependence (3.1) was derived for the case of equal transverse momenta. This condition is impractical and not 
the way the data were taken, so we extend the calculations to the case of unequal transverse masses. The parameters 
of the model are the critical temperature, T c (~ m^), and the space-like freeze-out hypersurface, t 2 — x 2 = t 2 = const. 
We assume a Gaussian transverse distribution exp(— r 2 ^/ R 2 ^) of hot matter in a pipe with effective area S± — TtR\. 



(rrii cosh#i, rrii sinhfj, pi), where pi is the transverse 
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The particles are described by their momenta fcf = (fc?, kj,p. 

momentum, 9i the particle rapidity in the beam direction, and mf = m 2 . + p 2 is the transverse mass. Further, let 
2a = 9\ — 6 2 , 29 = 8± + 9 2 and = p\__v 2 . Approximating the one-particle distribution by a Boltzmann form, the 
saddle point method yields from Eq. ( 2.11 ) (see Ref. Q): 



dN 
d9idpi 



TS±mi 



2irT c 
mi 



-mi/Tc 



(3.2) 



and 



: A formal change of variables can already be made at the operator level, Eq. (^6), but the calculation of the two-particle 
spectrum then becomes ambiguous 2 The problem is even more acute for photon emission: unphysical longitudinally 

polarized states are introduced 
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C 2 (fci,fc 2 )-1 



dN 2 



dB\dp\dB 2&V2 



dNi dNi 



dOidpi d92dp2 



- 1 



(3.3) 



7 f(Wx\Rx) rL ;, A , ,, . , exp<; - — 



[/l(z)/l(l/z)] 3 / 2 



x cos 



JL 

T, 



/i(z) sm — h ti( -) sm 



H{z) , ,1. F(l/z) 1 
/i(z) cos — ^ + /i(-) cos — ^— ^ -2 



3 

+ 4 



ff(z) + J ff(- 



+ G(z) + G(in , 



where fi = (mim^) 1 / 2 and z = (TO1/TO2) 1 / 2 , and we have introduced the functions 

1 



h[z) = 
tan H(z) = 



z z -F z {z- - f + 4F 2 sinh 2 a 

z 



2 1/4 

2/_2 „„„uo„.\2 



4F^(z^-cosh2o!)' 



(z 2 + cosh2a) 2 + F 2 (z 2 - X,) 2 

z z 

2F(cosh2a - z 2 ) 

;2_ jF 2( z _ l)2 +4j p2 s i n h 2 Q,' 



-,1/2 



f(\q ± \Rs_)=e-^ 2 



tan G(z) 



F(z 2 



1/z 2 



cosh 2a 



(3.4) 



(3.5) 



with F = tT c . Equation (3.2) is free of the limitation m\ — m2, and should be used directly to fit the model 
parameters, T c and r, without an intermediate Gaussian, and at various rapidities. The T c may then in turn be 
confronted with the slope of the one-particle spectrum. If, for instance, we find a dependence r(y) then it is natural 
to consi der the effect of a finite initial size via the 1-d Landau model as a next step. The visible size will be given 
by Eq. (3.1), with t — > r(y) [Q. Now, the relevant parameters of the model are the initial size and temperature, 
and the shape of the freeze-out surface. The general idea is that since there is no unique prescription to decode the 
interferometric data, the most effective approach is a "step-by-step" strategy - the complexity of the models used to 
fit the data should increase gradually in directions which are suggested by a thorough analysis of the simpler model. 

It has been conjectured H that kaons may escape from the source earlier and over a longer period of time than 
pions. Such a "dynamical" decoupling will influence the interferometric data. Suppose first that we have a sharp and 
common freeze-out of all particles; taking T c — 130 MeV and t = 30 fm, we obtain the curves marked "s" in Fig. [I]. 



As expected, we find that 7 



eft I 



1/2 • 



is independent of mi, and 7^ = jk- Suppose on the other hand that the 



kaons decouple dynamically from the pion fluid over some interval from tq to t; if tq = 7 fm (tT 3 = const then implies 
To = 180 MeV), and assuming for simplicity an unweighted average of emission function (2.11), (J(t)) t , we obtain 
the curves marked "g." Clearly, mj_ -scaling is violated for this type of emission. NA44 data |6|] gives no indication of 
such scaling violation, suggesting that the kaon fluid remains coupled to the pions right until their common (sharp) 
freeze-out. A systematic and parallel study of pion and kaon interferometric source sizes as a function of mi will 
provide a sensitive test for the temporal interval of freeze-out in hadronic matter at RHIC. 




0.4 



Ak||/mi /2 (fm- 1/2 ) 



FIG. 1. The correlation function C2, as a function of the rescaled longitudinal momentum difference. 



The dependence (3.1) may well apply at the SPS energies, since most of the matter will cool and reach the freeze-out 
stage due to the fast longitudinal expansion rather than the much slower transverse expansion, but at the AGS it is 
well known that transverse flow plays an important role. We illustrate the effects one might encounter by way of two 
unphysical, but exactly calculable examples: 
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2. Scale— invariant 2 d flow and the explosion of a long filament: To find a possible signature of transverse flow, 
consider the extreme case of a transversely expanding filament of length L R±. A convenient parameterization of 
the coordinates is x^ = (r cosh/3, r sinh/3 cos?/>, r sinh/3 smip, z), where (3 is the radial rapidity of the fluid cell. 
The temperature and velocity field may be written as r 2 T c 3 = const, and = (cosh/3, sinh/3 cos'0, 
sinh/3 sinV>, 0), while dE M = 
m z i = m z2 , where m z = ^pp\ 



^ t sinh/3 dzd/3dip. The emission function J(kx,k 2 ) may be calculated exactly for 
m 2 is the longitudinal mass. We obtain 



k 



dNi 
dk 

C 2 (k u k 2 )-1 



J(k,k) = 2itLt 2 T c 



(m z /T c ) 



Re 



-2m z {H-l)/T c 
ff6 



m z /T c 



m z 



m z H/T c + 1 



Q 2 



8m 2 



(3.6) 
(3.7) 



where Q 



-(h 



k 2 )\ 



H = 



1 



rpl 2 Q 2 

J-„T - 



T c tQ 2 



-,1/2 



2m? 



(3.8) 



Similar to the case of 1-d boost invariant flow, we obtain a plateau for cylindrical boost invariant expansion, but now 
in the radial rapidity distribution. The dependence of the spectral density on m z is noteworthy: the localization of 
the spectrum due to radial flow is more pronounced for greater m z . Particles with m z /T c 3> 1 spend almost all their 
thermal energy for the longitudinal motion and thus are strongly frozen into collective radial flow. For this reason 
we may expect the width of the correlator in the transverse direction to be defined by m z , rather than m±. This 
simple and well-understood dependence leads us to conclude that nit scaling of the transverse radius cannot be due 
to hydrodynamic expansion, and it is even doubtful that it is at all possible. 

In the case of small differences in radial particle rapidity, the exact result ( |3.7| ) may be simplified to obtain the 
same "radii" for the sideways and outward directions: 



Mc 1 

Rexp — ^"w ) Rcos — Ttt R-exp 

V m z l r T 



(3.9) 



where R eX p dominates the shape of the correlator for T c t > 1, while R cos dominates for T c t < 1. If the radial flow 
takes place against the background of a strong longitudinal expansion, which is likely in a more realistic scenario of 
the heavy ion collision, the distinctive m 2 -dependence in the data could be only washed out, but not replaced by an 
m^-dependence. 

3. Scale— invariant 3— d flow and the explosion of a point-like source: For the case of spherical expan- 



sion, — t (cosh/3, sinh/3 sin 9 cos ip, sinh/3 sin 6 sin?/;, sinh/3 cos 9), with tT c = const, 



x^/t, and 



dW = u^t 3 sinh z /3 sin6» df3d9dip. Once again, the emission function J(ki,k 2 ) may be evaluated exactly, and 
we obtain: 



dk 

C 2 (k u k 2 )-1 



4nr s T c K 2 ( Y ) , 



(3.10) 



Re 



K 2 {H) 



H 2 K 2 (m/T c ) 



1 



8m 2 



where H 



1 + Th 



Q 2 



T c tQ 2 

\ 

2m 2 



and K 2 is the modified Bessel function. We see that spherically symmetric boost-invariant expansion washes out any 
inhomogeneity of the local thermal spectrum. This is a distinctive feature of the spherical model. Approximating the 
correlator in the case of small differences in radial particle rapidity, we obtain 



C 2 {Q)-l 



1 



2 /4m 2 



l + 5T 2 r 2 Q 2 /2r 



Q\5T c t 

ros 1 — 2 (^T~ 

m z 2 



rriT . 



(3.11) 



All comments related to Eq. (3.9) apply here, too. Here, however, all three directions are physically equivalent and 
the mass of the particle becomes the only parameter that may re-scale the correlation radius. 



IV. CONCLUSION 



We have emphasized some important physical aspects of the theory of interferometry, and tried to isolate some 
controversial points. Strictly speaking, interferometry does not permit the initial data to be given in terms of semi- 
classical distributions, unless this description is augmented by a clear indication of the length scale that defines the 
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quantum states of the particles. Our conclusion is that the traditional operator approach, based on the precise 
definition of the particle states, provides a firm footing for the calculation of the two-particle spectra. 

Recent preliminary results by NA35 and NA44 seem to confirm the main predictions of interferometry for a scale- 
invariant 1-d flow scenario. They provide a multidimensional test of the source, and the fact that so many parameters 
coincide can hardly be accidental. Thus, we have strong evidence that at SPS energies a hydrodynamic regime 
develops, and that freeze-out takes place during a short interval. This collective behavior can be expected to occur 
at RHIC energies, and thus it is highly desirable to continue developing the formalism. 

We have extended previous calculations for a longitudinally expanding system to the case of unequal transverse 
momenta, and have shown that the HBT correlator may carry information about the time of the kaon emission. 
We further considered transversally expanding sources, and demonstrated how to change the parameterization of the 
correlator. Unfortunately, the attractive feature that the spectrum is localized in the boost-invariant solution for 
1-d flow is absent in 2- and 3-dimcnsional flows, even if the radial motion is very strong. This happens because the 
angular dependence of the Cartesian velocity components is weak, and practically means that one should perform 
all integrations exactly. While even moderate radial flow does obscure the true transverse source size, the spectrum 
is not sufficiently localized to allow one to obtain a simple formula for the correlator (at least in terms of standard 
variables). Numerical calculations are, of course, possible, but then it becomes difficult to "recognize" the model. 

We emphasize that while HBT for hydrodynamic sources is well understood as a physical phenomenon, choosing 
an adequate model to fit the data always poses difficulties. We impose the condition that the model should allow one 
to recognize it via a qualitative analysis. Only then can one hope to understand what parameters are responsible for 
the behavior of the correlator, and find their values by fitting the data. Practically, this requirement means that we 
must begin with an analytic expression for a solution of the relativistic hydrodynamic equations. This guarantees the 
consistency between the shape of the freeze-out surface and the velocity field. Unfortunately, analytic solutions for 
the case of three-dimensional expansion are not yet known. A reasonable analytic approximation will do, but to our 
knowledge, a suitable expression has not yet been derived. An approximate formula describing a realistic, expanding 
system at the freeze-out stage is an important problem for boson interferometry at RHIC. 
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